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Abstract 

The problems of escape from metastable state in randomly flipping potential and of 
diffusion in fast fluctuating periodic potentials are considered. For the overdamped 
Brownian particle moving in a piecewise linear dichotomously fluctuating metastable 
potential we obtain the mean first-passage time (MFPT) as a function of the po- 
tential parameters, the noise intensity and the mean rate of switchings of the di- 
chotomous noise. We find noise enhanced stability (NES) phenomenon in the system 
investigated and the parameter region of the fluctuating potential where the effect 
can be observed. For the diffusion of the overdamped Brownian particle in a fast 
fluctuating symmetric periodic potential we obtain that the effective diffusion co- 
efficient depends on the mean first-passage time, as discovered for fixed periodic 
potential. The effective diffusion coefficients for sawtooth, sinusoidal and piecewise 
parabolic potentials are calculated in closed analytical form. 
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1 Mean first-passage times in dichotomously fluctuating metastable 
potential 



The problem of escape from metastable state in fluctuating potential is of great 
importance to many natural systems, ranging from physical and chemical sys- 
tems to biological complex systems. The characteristic feature of all these 
complex nonequilibrium systems is that they are open systems, which are in 
contact with a fluctuating environment. A typical problem is the enhancement 
of stability of metastable states in fluctuating potentials due to the external 
noise [1,2,3,4,5]. The noise enhanced stability (NES) phenomenon was ob- 
served experimentally and numerically in various physical systems. This effect 
implies that a system remains in the metastable state for a longer time then 
in deterministic case and mean first-passage time (MFPT) has a maximum at 
some noise intensity. Specifically NES phenomenon was found in the transient 
dynamics of a periodically driven overdamped particle in a noisy cubic poten- 
tial [4], in a tunnel diode by investigating the escape time in a periodically 
driven metastable system in a strong forcing regime [1]. More recently it was 
found that several properties of the noise-induced stability in one-dimensional 
maps and in fluctuating bistable potentials are related to NES effect [6]. The 
noise induced slowing down exhibited by the mobility of an overdamped par- 
ticle moving in a periodic potential and in an inhomogeneous medium [7], the 
noise induced order in one-dimensional map of the Belousov-Zhabotinsky re- 
action [8], and the one-peak structure in the MFPT of a bistable kinetic model 
driven by correlated noises [9] are akin to the NES phenomenon. Previous the- 
oretical papers analyzed NES phenomenon in systems with metastable and 
unstable states of fixed potential or periodically driving potential [1,2,3,4,6]. 
However the model of randomly switching metastable state is more realistic 
when we describe the generation process of the carrier traps in semiconductors. 
Let us consider one- dimensional overdamped Brownian motion in randomly 
switching potential profile 



where x{t) is the displacement in time t, ^{t) is the white Gaussian noise 
with zero mean and {^{t)^{t + t)) = 2D6{t), ri{t) is Markovian dichotomous 
process which takes the values ±1 with the mean rate of switchings u. We 
assume that potential U{x) has reflecting boundary at the point x = and 
absorbing boundary at the point x = b {b > 0). Exact results for MFPTs 
for Brownian diffusion in switching potentials were first derived in [10]. Here 
we use equations from [10] but different conditions at the reflecting boundary 
[11]. Thus following ref.[10] we obtain the coupled differential equations for 
the MFPTs of our system (1) 
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Fig. 1. The piecewise linear potential U{x) with metastable and unstable states. 

DTI + [a-u' (x)] r; + 1/ (T_ - r+) = -i, 

DT" -[a + U' (x)] T'_ + v (T+ - T_) = -1. (2) 

Here T+(x) and T_(x) are respectively the MFPTs for positive 77(0) = +1 
and negative 77(0) = —1 initial value of dichotomous noise. In accordance with 
ref. [11] the conditions at reflecting boundary x = and absorbing boundary 
X = b are (0) = and T± (b) = 0. Let us introduce for convenience two 
new auxiliary functions T = (T+ + r_)/2 and 6 = (T+ — TS)/2. So that from 
Eqs. (2) we obtain 



DT" -U' (x) T' + ae' = -1, 

DO" - U' (x) 6' + aT' - 2^9 = 0. (3) 

with the following boundary conditions: T' (0) = 9' (0) = and T (b) = 9 (b) = 
0. By eliminating T(x) from Eqs. (3) we arrive at third-order linear differential 
equation for 9{x) 



a'" (•^) n" 



(x) U" (X) 



D 



9' 



2uU' (x) 



a 



(4) 



where 7 = a? / D"^ + 2z//D. Let us assume that in Eq. (1) U{x) = k{L — x) ■ 
l(x — L), where l(x) is step function and < L < 6. As we can see from 
Fig. 1 we have a metastable state for 77 (t) = —1 and an unstable state for 
r/(t) = +1. We will consider the MFPT T+(0), i. e. MFPT for initial unstable 
state with starting position of Brownian particle at the origin. In the absence 
of switchings and thermal diffusion the dynamical MFPT T+(0) equals 



We can solve Eqs. (3) and (4) separately for regions < x < L and L < x < b. 
After some algebra we obtain finally the MFPT T_|.(0) as 



T+ (0) =T{0)+e{0) = c,^ + c,- -4frr (6) 

The unknown constants Ci and C2 can be determined from the boundary con- 
ditions and the continuity conditions of the functions 9{x), 6'{x), T{x), T'{x) 
at the point x = L. 



2 Analysis of noise enhanced stability phenomenon conditions 



By expanding the Eq. (6) in power series on D, for very low intensity we 
obtain 



2L vL'^ b-L N D , ^ , 

<°) T + 1^ + — - {l-e-)+-,F(sJ},.) (7) 



where 



^ ^ (1 + /3) (1 - + 2(1-/52) V 



2(l+/3) (3) 13^ 

and s and u are dimensionless parameters 



The sign of the function -F(s, /3, cu) in Eq. (7) determines the condition to ob- 
serve the NES phenomenon in considered system. We can write this condition 
by the following inequality 



F(5,/?,a;)>0. (9) 

Let us analyze the structure of NES phenomenon region on the plane u;). 
In the case of very slow switching z/— >0(ci;— >-0, s— >0) Eq. (9) takes the 
form 
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Fig. 2. The shaded areas are the region of the plane {uj,P) where the NES effect 
takes place. 



^> 0,802; uj< , 

6 (1 - /52)2 - 2/5 (1 - /?2) + /52 (3/52 - 1) - 1) 

In the case of /3 ~ 1 we obtain from Eqs. (8), (9) and (10) 



(10) 



"<67^' 2 + 2(^-^)<"<2(r3;^- 

In Fig. 2 are shown two NES regions (shaded areas) on the plane (/5, u). As we 
can see from Fig. 2, the NES effect occurs only at the values of /5 near 1, i.e. 
at very small steepness k — a = k{l — (3) of the reverse potential barrier for the 
metastable state (see Fig. 1). This NES region is different from that obtained 
for the case of periodical dichotomous force [2]. On our opinion, this narrowing 
of NES phenomenon area in comparison with the case of periodically driven 
metastable states is due to the uncertainty of metastable state starting time. 
The NES effect disappears when we choose the absorbing boundary at the 
point y = L. 



3 Acceleration of diffusion in fluctuating periodic potentials 



The problem of diffusion in a periodic potential is related to a rich variety 
of physical situations such as diffusion of atoms in crystals, synchronization 
of oscillations, fluctuations of a Josephson supercurrent and so on. An exact 
expression for diffusion coefficient was independently obtained in the over- 
damped limit by various methods for an arbitrary periodic potential and for 
sinusoidal potential [12]. Recently the mean velocity and effective diffusion co- 
efficient of Brownian particle moving in a tilted periodic potential have been 
found in [13]. The case of diffusion in fast fluctuating periodic field was first 
investigated in ref. [14], where the exact result for diffusion coefficient in the 
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sawtooth potential was derived. In this section we generahze this result to 
the case of arbitrary potential profiles (see also [15]). We will find that the 
problem of calculating an effective diffusion coefficient in fluctuating periodic 
potential reduces to mean first-passage problem similar to the case of fixed 
potential profile [16]. We consider an overdamped Brownian particle in fluc- 
tuating periodic potential U (x) whose dynamics is governed by the Langevin 
equation 



| = -^.cw+ew, (11) 

at ax 

where x{t) is the displacement in time t, ^ (t) and ( (t) are statistically inde- 
pendent Gaussian white noises with zero means and intensities 2D and 2Z)^ 
respectively. Further we assume that potential U {x) is even function with pe- 
riod L and place the origin in one of the potential minima. Following ref. [12] 
we determine the effective diffusion coefficient as the limit 



D^ff = hm (12) 



because we have {x{t)) = 0. From Eq. (11) the Fokker-Planck equation for 
probability density W {x, t) ^ is 



dW d , d , d'^W 

We can choose arbitrary initial condition for Eq. (13), since we look at the 
asymptotic behavior of mean square coordinate of Brownian particles. We 
place therefore all Brownian particles in the origin a.t t = 0: W [x, 0) = S (x). 
Because of periodicity of the potential, the diffusion process can be coarsely 
conceived as consecutive transitions of Brownian particle from points of poten- 
tial minima Xm = rnL to nearest neighboring points Xm±i- The transition time 
represents the escape time over left or right absorbing boundaries x = Xm±i 
for particle starting from the point x = Xm, i-e. the random first-passage time. 
Thus, we can consider as in [13], the discrete process 



n{0,i) 

m = E (14) 
k=0 

where are random increments of jumps with values ±L and n(0, t) de- 
notes the total number of jumps in the time interval (0,t). In the asymptotic 



^ We interpret the stochastic differential equation (11) in Stratonovich's sense. 
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limit t oo the "fine structure" of a diffusion is unimportant, and the ran- 
dom processes x (t) and x{t) become statistically equivalent. In particular, 
{x^ (t)) ~ {x^{t)). Further, the random increments Qi and the random times 
Tj between jumps are statistically independent of one another because of the 
markovianity of x (t) and have the same probability distributions W (q) and 
w (r) respectively. Because of the symmetry of potential U [x) the probability 
density W (g) reads 



W{q) = ^[5{q-L) + 5iq + L)]. (15) 
Let us calculate the second moment {x^{t)) from Eqs. (14), (15) 



i\t)) = Y.Pn it) EE M = Y.Pn it) E {qI) = (n(0, t)) , (16) 
n=0 k=Ol=0 n=0 k=0 

where P„ (t) is the probability of n jumps in the time interval (0,t). On the 
other hand in the limit t — > oo we have that 



(n(0,t))^-, (17) 

T 

where r = (tj) is the MFPT for Brownian particle with initial position x = 
and absorbing boundaries at a; = ±L. After substitution of Eqs. (16), (17) in 
Eq. (12) we arrive finally at the exact result 



Deff = ^, (18) 

which was first derived in ref. [16] for fixed periodic potential. In fluctuat- 
ing periodic potentials therefore the calculation of diffusion coefficient -De// 
reduces to mean first-passage time problem. According to Eq. (13) we must 
solve the following equation with conjugated kinetic operator 



Dt" (x) + DrU' (x) 4- \U' {x) t' {x)] = -1 (19) 
dx 

and boundary conditions r {~L) = 0, r (L) = 0. Then we must put x = 0, 
i.e. find r = r (0). Solving Eq. (19) with these boundary conditions, we finally 
obtain the following exact formula for effective diffusion coefficient of Brownian 
particle in fast fluctuating periodic potential U (x) 
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D 



dx 



D^[U'{x)Y/D 



(20) 



As evident from Eq. (20), -De// > D for an arbitrary potential profile f/(x), 
therefore diffusion of Brownian particles accelerates in comparison with the 
case U {x) = 0. This result fully confirms the assumption previously proposed 
in [14]. We emphasize that the value of diffusion constant does not depend on 
the height of potential barriers, as for fixed potential [12], but it depends on 
its gradient U' (x) . We can explain the phenomenon of diffusion acceleration 
directly from Eq. (18). The potential barriers change places through a random 
modulation and Brownian particles move from point x = to point x = L 
more rapidly in comparison with free diffusion case, i. e. in the average particles 
move downhill for the most part of the distance. Let us consider particular 
shapes of potential U (x) . 

(a) For the sawtooth profile U (x) = 2E \x\ /L at |x| < L/2 we easily obtain 
the Malakhov's exact result [14] 



D,ff = D + Dt 



4^2 
L2" 



(21) 



(b) For sinusoidal potential U (x) = Dsin^ (ttx/L) from Eq. (20) we have 



D, 



7r2D(l + 72) 



4K2 




(22) 



where K (k) is the complete elliptic integral of the first kind (0 < A; < 1). We 
derive from Eq. (22) at small strength D^ of modulating white noise (7 <^ 1) 
the following approximating expression 



D^fj ^ D + Dc^. 

This formula coincides with approximate result [14] obtained on the coarse 
assumption of Gaussian probability density W (x, t) although the real prob- 
ability density is multi-modal. In opposite case 7 1 using the asymptotic 
formula for elliptic integral 



we find from Eq. (22) 
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Deff ^ ^ ^' 23 

411127 In^D^ ^ ^ 

According to Eq. (23) the effective diffusion coefficient increases with intensity 
Dc^ of modulating noise but more slowly than linear law (21). 
(c) Finally for piecewise parabohc periodic potential profile 



Uix] 



E 



8E{x/L)\ \x\<L/4: 
- 8 {x/L - l/2f] , L/4:<x < 3L/4 



we get from the exact formula (20) the following result 



^e//= , . ^=—\\^- (24) 



In^ [m + vTr+"?r?^ 



At comparatively small intensity D^ (m <^ 1) we obtain 



that is quite similar to the formula obtained for sinusoidal potential. Moreover, 
the dependence of effective diffusion coefficient on large D^ is similar to the 
law (23) for sinusoidal potential profile. 
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